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^ Abstract 

Using holography, we study the entanglement entropy of strongly coupled field theories perturbed by 
^— H operators that trigger an RG flow from a conformal field theory in the ultraviolet (UV) to a new theory 

in the infrared (IR). The holographic duals of such fiows involve a geometry that has the UV and IR 
regions separated by a transitional structure in the form of a domain wall. We address the question 
I of how the geometric approach to computing the entanglement entropy organizes the field theory data, 

I exposing key features as the change in degrees of freedom across the fiow, how the domain wall acts as a 

• • UV region for the IR theory, and a new area law controlled by the domain wall. Using a simple but robust 

_ 1^ model we uncover this organization, and expect much of it to persist in a wide range of holographic RG 

flow examples. We test our formulae in two known examples of RG flow in 3+1 and 2+1 dimensions 
l] that connect non-trivial fixed points. 
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1 Introduction 



A useful probe of the properties of various field theories that has received increased interest in recent times 
is the entanglement entropy, with applications being pursued in diverse areas such as condensed matter 
physics, quantum information, and quantum gravity. One of the main motivators, in the context of strongly 
coupled field theories (perhaps modeling novel new phases of matter) , is that the entanglement entropy may 
well act as a diagnostic of important phenomena such as phase transitions, in cases where traditional order 
parameters may not be available. 

Within a system of interest, consider a region or subsystem and call it A, with the remaining part 
of the system denoted by B. A definition of the entanglement entropy of A with B is given by: 



where is the reduced density matrix of A given by tracing over the degrees of freedom of B, = Tre(p), 
where p is the density matrix of the system. When the system is in a pure state, i.e., p — the 
entanglement entropy is a measure of the entanglement between the degrees of freedom in A with those in B. 

It is of interest to find ways of computing the entanglement entropy in various strongly coupled 
systems, in diverse dimensions, and under a variety of perturbations, such as the switching on of external 
fields, or deformations by relevant operators. A powerful tool for studying such strongly coupled situations 
is gauge/gravity duality, which emerged from studies in string theory and M-theory. The best understood 
examples are the conjectured AdS/CFT correspondence and its numerous deformations [J [51 [31 13] (See 
e.g., ref.[S] for an early, but still very useful, review.) There has been a great deal of activity for over a 
decade now, applying these tools to strongly coupled situations of potential interest in condensed matter 
and nuclear physics, for example. Fortunately, there has been an elegant proposal^ [7] for how to compute 
the entanglement entropy in systems with an Einstein gravity dual (or, more generally, a string or M-theory 
dual in the large N limit and large t' Hooft limit), which provides a new way to calculate the entanglement 
entropy using geometrical techniques (for a review see ref.[8j. In an asymptotically Anti-de Sitter (AdS) 
geometry, consider a slice at constant AdS radial coordinate z = a. Recall that this defines the dual field 
theory (with one dimension fewer) as essentially residing on that slice in the presence of a UV cutoff set by 
the position of the slice. Sending the slice to the AdS boundary at infinity removes the cutoff (see ref. [2| for 
a review). On our z = a slice, consider a region A. Now find the minimal-area surface 7,4 bounded by the 
perimeter of A and that extends into the bulk of the geometry. (Figure [T] shows examples of the arrangement 
we will consider in this paper.) Then the entanglement entropy of region A with B is given by: 



Tr^ {PA Inp^) 



(1) 



Sa 



Area(7^) 
4Gn 



(2) 



where Gn is Newton's constant in the dual gravity theory. 
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(a) The strip. (b) The disc. 



Figure 1: Diagrams of the two shapes we will consider for region A. This is the case of AdS4, and here, z denotes 
the radial direction in AdS4. In one dimension higher we will generalize these shapes to a box and a round ball, and 
in one dimension fewer, we will consider an interval. 



This prescription for the entropy coincides nicely with various low dimensional computations of the 
entanglement entropy, and has a natural generalization to higher dimensional theories. Note that there is 
no formal derivation of the prescription. Steps have been made, such as in refs. [HI [in]i but they are not 
complete. However, there is a lot of evidence for the proposal. See e.g., refs.[TTl[T21[I31[Ill[ini[IS]- A review 
of several of the issues can be found in ref . 8J . Further progress has been made recently in ref. |17j . 

In this paper we shall assume that this holographic prescription does give the correct result for the 
entanglement entropy in systems with gravity duals, and proceed to examine the interesting question of how 
the entanglement entropy behaves when a system is perturbed by an operator that triggers a Renormalization 
Group (RG) flow. For simplicity, we will work with flows that connects two conformal field theories, and we 
will consider (for concreteness) a four dimensional example and one in three dimensions. Such examples are 
extremely natural to study using holographic duality since (at large N) it is possible to find geometries that 
represent the full flow from the maximally supersymmetric theory to theories with fewer super symmetries. 
(This was first proposed in refs. [lH [19], and several examples have since been found.) Flow between field 
theory fixed points correspond to flows between fixed points of the supergravity scalar potential. The 
examples we will study begin with the four dimensional case of the flow[5ni HH [22] to the Leigh-Strassler 
point [23l [21] , which results from giving a mass to one of the N — 1 chiral multiplets that make up the M — A 
Yang Mills gauge multiplet. We then continue with the three dimensional generalization of it discussed in 
ref. |2S]. The gravity dual of the four dimensional flow connects AdSs x at the r = +oo extreme of a 
radial coordinate r to AdS^ x Ai^ at r = — oo, where the space Ai^ results from squashing the along the 
flow. There are two of the 42 supergravity scalars switched on at the latter endpoint, and correspondingly 
the characteristic radius of the AdSs in the IR is larger than that of the UV theory: The gravity dual for 
the three dimensional flow has related structures, this time connecting an AdS4 x 5*^ UV geometry to an 
AdS4 X A^7 in the IR, where Mj results from squashing the S'^ along the flow. 
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Before studying the specific examples, however, we step back and try to anticipate some of the 
key physics that we should expect from the entanglement entropy in this type of situation, more generally. 
Generically, holographic RG flow involves a flow from one dual geometry in the UV to another in the IR, 
separated by an interpolating region that can be thought of as a domain wall separating the two regions. 
The key to understanding the behaviour of the content of the holographic entanglement entropy formula is 
to then understand how the computation incorporates the structure of the domain wall, and how the field 
theory quantities it extracts are encoded. To anticipate how to mine this information, we do an analytic 
computation of the proposed entanglement entropy ^ in an idealized geometry given by a sharp domain 
wall separating two AdS regions with different values for the cosmological constant. Working in various 
dimensions (AdSs, AdS4, and AdSa, pertaining to fiows in four, three, and two dimensional field theories), 
we find a fascinating and satisfying structure, seeing how the entanglement entropy tracks the change in 
degrees of freedom under the flow, and several other features. We expect that these features will be present 
in a wide range of examples, and we confirm our results in the examples mentioned above. 

The outline of this paper is as follows. In section 2 we carry out the study of the entanglement entropy 
in the presence of the idealized {i.e., sharp domain wall) holographic RG flow model, and discover how the 
physics is organized in the results. Then, ready to study examples, we review the four dimensional Leigh- 
Strassler RG flow of interest, and its dual AdSs flow geometry in section 3. We explicitly solve (numerically) 
the non-linear equations that define the geometry and scalars in the interpolating dual supergravity fiow. 
We then compute the entanglement entropy and extract the physics, comparing to our predictions from 
section 2. Section 4 presents the analogous studies for the three dimensional field theory, with the AdS4 dual 
fiow geometry. We end with a discussion in section 5. 

2 Entanglement Entropy and a Sharp Domain Wall Model 

As mentioned in the introduction, the generic holographic RG flow involves a flow from one dual geometry in 
the UV to another in the IR, separated by an interpotating domain wall. In all examples, understanding the 
behaviour of holographic entanglement entropy, as proposed in equation ([2]), requires us to understand how 
the area formula incorporates the structure of the domain wall in terms of field theory quantities. So we start 
by doing an analytic computation in an idealized geometry given by a sharp domain wall in AdS. In general, 
the location of the wall, and its thickness, are determined by field theory parameters corresponding to the 
details of the relevant operator - for example, in the case of the Leigh-Strassler flow and its generalization we 
later study, the detail in question is the bare value of the mass given to the chiral multiplet. A sharp domain 
wall is of course not a supergravity solution, and falls somewhat outside the usual supergravity duality to 
any (large N) theory, but nevertheless is a clean place to start to capture how the physics is organized. We 
expect it to capture a great deal of the key physics of holographic RG flow, as regards how the entanglement 



4 



entropy formula works. 

We use the following background metric: 

ds^ = e2^('-) {-dt^ + dx^) + dr^ , (3) 

with 

Here — oo < r < +00, and x is either four, three, or two coordinates (the spatial coordinates of the dual 
field theory), depending upon whether we are in AdSs, AdS4, or AdSa, the cases we will consider. Also, 
i?iR > i?uv- The length scale of AdS on either aide of the wall is set by i?uv in the UV at r > row and 
RiR in the IR at r < Rdw- 

2.1 The Ball and AdSg. 

We begin by studying a region A in the three spatial dimensions which is a round ball of radius £. Using a 
radial coordinate p in the spatial dimensions, the area of the surface, 7, that extends into the bulk is given 
by: 

Area = An j' dp p'e^^^''^ (l + e-^^^^\'{pf) , (5) 

where the function r{p) defines the onbcdding. We can calculate the equations of motion that result from 
minimizing this "action," and we find that the solution is given by: 



r{p) 



2 - e^^fe^) ' P^P'^^ 



-^In 

V -"uv J ,„s 



where row is the position of the domain wall in the AdS radial direction, and pow is the spatial radial 
position where r{pDw) = Tdw and is given by: 

pIw = + - ii?rve"'^ • (7) 

Note that, rather than integrating out to r = +00, we integrate out to large positive radius ruv> defining 
our UV cutoff, with small e defined by: 

ruv = --Ruvln( 1^ ) • (8) 



Note that with the solution given for we are assuming that i is larger than a critical radius icr such 
that our surface extends past the doman wall into the second AdS region. The critical radius ^cr is given by 
setting pdw = in the above: 

C = ii^ve"'^ - • (9) 
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Substituting the solution back into equation ([5|, we can analytically calculate the area of our minimal 
surface, 7yi, and hence the entanglement entropy via equation ([2]). This gives a long expression that we will 
not display here. For our purposes it is enough to first expand the area for small e: 

Area 



7?3 



47r 



^2 



In 



-«UV 



[l-21n(2)] 



1 



Rl^e "«uv + -i?3jytanh-i 



P - Rlye~^^ 



P + Rf^e- 



?3 
2 



^tanh-i 



\ £2 + ^: 



2 -2- 



We find it useful to rewite it in a suggestive way: 



Area 
47r 



^2 , /£ 
^ +ln 



%[l-21n(2)] 



tanh 



0(6) 



_9 LDHl 



- £2 



(10) 



where 



-RlR 



2 
2 



^2 - £2^ + Rl^er^-i^ 



tanh-1 



\ P^P„ + Rl^e~^^ 



R 



0(e) 



(11) 



(12) 



There are a number of notable features of this expression. First, we see the results from pure AdSs in the 
first line. There, we see the usual UV divergent terms and the ^-independent constant that results from the 
fact that the ball preserves some of the conformal invariance of AdSs. Second, the terms that have as 
coefficients (the last two lines) always have t"^ appearing in the combination: 



£2 



(13) 



We are tempted to interpret this £ as the efi'ective ball radius as seen in the IR, as opposed to the simple i 
seen in the UV. Furthermore, the combination: 



-RiRe 



(14) 



appears in a manner analogous to how the UV cut-off e appears. (This might not be clear in our e expansion 
of the above equation. One way to see that it does appear as e does is to look at equation ([6])). Now e is not 
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necessarily small, but we will see that it is useful to think of it as the cut-off in the IR theory. With these 



observations in mind, we rewrite the last three lines of equation (|11|) as follows: 



7?3 



2 



i?uv tanh 




9 ?2 



Mil 
2 



tanh-i 




(15) 



Let us focus on the terms proportional to If we expand these terms assuming that e/t « 1, i.e., the 
effective length in the putative IR theory is larger than the IR cutoff, which also means that the length is 
such that the surface extends very far past the wall into the IR AdS space, we get: 



R 



IR 



[l-21n(2)] +0(g/^~) 



(16) 



Pleasingly, this is exactly the result we would have obtained if we were purely in the IR theory! 

So far therefore, we have seen how the entanglement entropy formula encodes key behaviours of 
both the UV and the IR theories, in terms of the appropriate scales, e/£ and e/l. The boundary of AdSs 
at r = +00 is the UV region and the quantities of the UV theory appear accordingly. From the point of 
view of the IR theory, the domain wall acts (for e/£ small) as the effective UV region, with e/l acting as the 
effective regulator. 



We are left with understanding the first two terms in equation (15). These two terms mix the 
properties of the UV and the IR regions, and are more subtle. We associate them with the region around the 
domain wall, which connects the UV and IR regions (through and abrupt change in our idealized example). 
It is prudent to try to understand the role of these terms toward the end of the flow, and so we do a large £ 
expansion of them, giving: 



2 'W 



1 



i?uv tanh 



In 



7?3 



+ ^[l + 21n(2)]+0(l/^) 



(17) 



So we see that these terms give contributions very analogous to our UV and IR results, where here the 
reference scale is played by i^r- (Note that the constant term is actually different than the UV and IR 
constant terms' form.) At fixed e or e, we may think of this as a new set of divergences. 

Now that we have an understanding of the contributions of the various pieces to the area, we combine 
everything together again and consider the large i (and small e) expansion: 



Area 






47r 


2 










2 












2 










2 





In 



^IR 



In 



-"-IR 



2 



£Hjv 

£2 



1 



4 



1 n3 "cr 
2^IR ?2 



2 

£2 

^2 
2 1 
2 



r!^H2) + o 



^In 
2 

1 



^In 



-"^IR 



-"-UV 



R^ 

^In 



^i??Rln(2)+0( ^, e 



(18) 
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The first key result here is that we no longer have a \n{£) scaling associated with the UV theory. The 
remaining ln(£) dependence has a coefficient that is only associated with the IR theory and that is independent 
of the domain wall. In a non-RG flow scenario, the coefficient of such a term is determined by the central 
charge of the theory (see e.g., refs.[7| US]), but here we see that the coefficient has shifted from its UV value 
(associated with the UV central charge) to its IR value (associated with the IR central charge). The second 
thing to note is that the area law associated with the UV cut-off (the ffi-st term) is joined by a second 
area law. Its coefficient is sourced by the details of the domain wall. For clarity, we display this term here: 



2 I f2 i2 ) - 2 



i?jp,e -i?uve «uv . (19) 



We expect this new area law to be a robust feature of RG flow geometries, but anticipate that the coefficient's 
precise form will be different as we move away from the thin wall limit we are in here. The above result 
predicts that the coefficient grows more positive as tdw is pushed to the UV. In realistic RG flows, while 
the domain wall position and sharpness cannot be varied arbitrarily, it is expected to get thinner toward the 
UV and so at least in that regime we should recover positivity. Finally, the constant terms in the last line 



of equation ( 18 1 are a mixture of both the UV, IR, and domain wall physics. 



2.2 The Disc and AdS4. 

We can repeat the same procedure for AdS4, pertaining to RG flows in 2 + 1 dimensional theories. As our 
system A we consider a circular disc of radius I. The solution for the surface embedding are exactly as in 
equation (|6|. We can calculate the minimal area and expand for small e to get: 



Area o £ o t>»' ^dw / „ „ o t>»' „„ o ^nw „, , , , 

— =Rl^-- R!n - i?uve £ + R,^e .J p - «uv + Rf^e^ + 0(e) . (20) 

We see the reappearance of many of the key players that we saw in the AdS^ case, such as £, £cr and e, 
appearing in similar types of term. For £ ~ £cr, we recover the pure UV result (proportional to £/e) and also 
the constant —R^y, the constant ensured by the fact that the disc preserves some conformal invariance, as 
expected. For large £, we have: 

^"^"^ Rl/- - RlyS- + RlJ, - + 0(e, II £) . (21) 



2n "'e ""U^r e 
So in the AdS4 case, the constant term shifts from its UV result to its IR result —Rjj^. Again, in addition 
to the usual UV area law (proportional to £/e), we have a new area law controlled by the domain wall: 



p2 n2 
^IR ^UV 



(22) 



which should be compared to the example from AdSs in equation (19). The same comments we made for 
the new area law there apply here: It is not necessarily positive, but we expect it to get more positive as the 
domain wall is sent to the UV, where generically it gets thinner. 
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2.3 The Case of AdSg. 

Next we consider the case of AdSa, pertaining to flows in 1 + f dimensions. We use a spatial interval of 
length 21 for our region A. The area is given by: 



Area 



^iT!uv In + i?uv ln(2) - i?uv tanh" 




2 ''PW' 

Rive '«uv 



+i?iR tanh ^ 



2 e '^uv +i?2j^e 



0(6) 



In the large £ limit, this gives 
Area 



-i?uv In 



uv 



i?iR In 



i?iR,ln(2) + 0(e,l/£) , 



(23) 



(24) 



where ^cr and e are defined in equations (12) and (14) respectively. So again we see that the universal 
coefRcient (in front of the natural logarithm) becomes the IR factor in the large £ limit. The IR cutoff 
replaces the UV cutoff just as observed before. 



2.4 The Strip and AdS4. 

We next consider an area A that is a strip in AdS4, to compare our results for the disc. We take the strip 
to be of finite width £ in the x direction, and of length L in the remaining direction, which will be taken to 
be large, making an infinite strip. The area is given by: 

/2 



Area = 2L 



(25) 



Since there is no explicit dependence on cc, there is a constant of motion in the dynamical problem associated 
to minimizing the area. However, we must be careful since the constant of motion on either side of the domain 
wall is not the same: 

°2^('-) I p2H^ ^ ^ 



2-1/^ 



^l+e-2A(r)^/(2.)2 



r < row 



(26) 



On the IR side, the constant is simply given by r'{x) = 0, which occurs at a radial position we will denote 
as r*. The constant on the UV side is determined by asking that r'[x) ~ when r.^, = row^ which is the 
critical situation before our embedding enters the IR AdS. We can in turn calculate the area and length in 
terms of and expand for small e: 

Area 



2L 



— - i?uve «uv 



e 

-e "IR 



1 



^ ^Ri^T (I) 1 

^TR — —i 



— e -"uv 
3r(|) 3 



1374: 

2'4'4'^ 



1 g ^DW 

3 



(27) 
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-36 ''^- +'"™ i?iR2Fi f «iH j. (28) 
Here, we see the appearance of the Gauss hypergeometric function: 

^' " - fm^) I (i-Jr ■ 

The large hmit corresponds to taking —00, which gives us: 

Area ^nv „ '•-p'^ „ '•pw ^/ , , ^ / x 

= -i?uve«uv +ii'iRe «iR +0(e,-l/n) . (30) 

So we see that that the constant term here is exactly the new area law's coefficient that we saw in the disc 



case, in equation (22). Again, far enough in the UV, for large enough mass, our analysis suggests that this 



coefficient is positive. 

2.5 The Box and AdSg. 

Returning to AdSs, we consider for region A a box in AdSs, in order to compare to the round ball we studied 
before. Here the finite width is again £ and the two other sides are of length L, which we again take to be 
large. The computation proceeds in a similar way. The area gives: 



Area R^y 1 A e '-^"""' , '^^xh e "* „''°'^ 

= - ^-^uve «uv Y 1 - e «uv + -Rme «m yi~e 

2^y/nRmr{l) 1 _4Lam_^_Q_J:^ /I 2 5 e -'-''ow 

8T{1) 8 V2 3' 3' 

1 _4Tpt+6^ /I 2 5 el±_fnw\ 

+ ge «m ^'' ^^2 ^2' 3' 3'^ j ' ^ ' 



and taking the — > —00 limit gives: 



Area _ i?^v 1 n '^iF^ , ^ u 2 



i?uve^uv + R^j^e^ Rm +0{e,-l/r^) . (32) 



2L2 2e2 2 2 

Again, we have that the constant term has the same coefficient as the new area law term for the ball case. 



as seen in equation (19). 



3 The Four Dimensional Holographic RG Flow 
3.1 The Holographic Dual Gravity Background 

In field theory terms, the RG flow is defined by an A/" = 1 supersymmetric deformation of the A/" = 4 
supersymmetric Yang Mills theory given by introducing a mass term for one of the chiral multiplets.This 
relevant deformation causes the A/" = 4 theory to flow to an A/" 1 fixed point in the IR called the Leigh- 
Strassler fixed point [321 1131 [Ml- For the SU{N) theory at large N, there is an holographic dual of this 
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physics [20]) represented by a flow between two five dimensional anti-dc Sitter (AdSs) fixed points of A/" = 8 
gauged supergravity in five dimensions. One point has the maximal SO(6) symmetry, and the other has 
SU{2) X [/(I), global symmetries of the dual field theories. The relevant five dimensional gauged supergravity 
action is [5D1 HH [H] : 

S = / d^xV^ (i? - 2 {dxf - 12 {daf - AV) , (33) 

with 

where the superpotential W is given by: 

W=^ (cosh(2x) {p' - 2) - (3p« + 2)) , (35) 

with p — exp(a). The scalar field x is dual to an operator of dimension three in the field theory while the 
scalar field a is dual to a dimension two operator: 

4 6 

a:^Tr(0,(/.,)-2^Tr((/.,0,) , x : Tr (A3A3 + [(^2, ^^s]) + h.c. , (36) 

i—l i— 5 

where tpk = 4>2k-i + i4>2k, k = 1, . . . , 3. Here 4>i (* = I7 • ■ • 7 6) are the six scalars in the A/" = 4 multiplet, 
and the (three of that adjoint multiplet's four fermions) are Af = I partners of the cpj., forming the three 
chiral multiplets. This combination of operators is exactly what is needed to reproduce the deformation. 
The geometry in five dimensions, of domain wall form, can be parametrised in the following manner: 

dsli = e^^'^''^ {-dt^ + dxl + dxl + dxl) + dr^ . (37) 

The supergravity equations of motion yield the following flow equations: 

da _ e°' dW _ 1 f e^" (cosh(2x) - 3) + cosh(2x) + 1 
dr 

dx _ I dW _ I ( (e^" - 2) sinh(2x) ' 
dr 




dA _ 2 1 cosh(2x) (e^" - 2) - (3e 



(38) 



dr 3R 6R e^" 

In these coordinates, the UV is at r — )• +00 and the IR is at r — > cx), as in earlier sections. In cither limit, 
the right hand side of the first two equations vanish, and the scalars run to specific values (a = 0, x = at 
one end, a — g ln2,x = ^ In 3 at the other), while A{r) becomes in the UV and ^-^-^ in the IR, defining 
an AdSs in each case, and hence a conformally invariant dual field theory at each end. This is a fat, smooth 
version of our simple thin domain wall model of the previous sections. Here R\jv — R and i?iR — 3i?/2^/^. 
To study the UV behavior of the fields, we find it convenient to define a coordinate z given by: 

S = e-^l^ , (39) 
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and we find the asymptotic behavior of the fields near z ^ (UV AdS boundary) : 
X{z) = s(^aQ + S^ (^-~al-4:aoai+\n{S)^al^ +0{S*) 
a{z) = f2 (^ai+ln(z) (^-^a^^^ +0(z4) , 

A{z) = -\n{2)+Ao-^al2^ + OiS''). (40) 
The constant oq is related to the mass of the $3 multiplet via [2S] : 

"^3 = • (41) 
To study the IR behavior of the fields, we define a coordinate u given by: 

u = e^'-/^ , (42) 

where A — -. The asymptotic (near u — 0) behavior in the IR is given by: 

X{u) = ^In(3) + fi6o + 0(u2) , 

a(u) = llni2) + u(^^^^^bo + 0{u^) , 

A{il) = —} ln(u) + Bo + 0{u^) . (43) 

V7 - 1 

3.2 Numerically Solving for the Flow 

To solve the fiow numerically (as we will need to do in order to compute the entanglement entropy), it is 
convenient to work with a coordinate: 

X = , (44) 

and employ a shooting method to solve the equations. To shoot from the IR we take Xmax — 10^, and 
towards the UV we take x = e, where we use e = 10^^^. To get good numerical stability, we define new 
fields {(3,ri), 

a{x) = xf3{x) , x{x) = x^^^-qix) , A{x) = ln(x) + a{x) , (45) 
such that near the AdS boundary, the leading behavior of these fields is given by: 

f]{x) = ao + 0{x) , f3{x) — ai — -a^ \n{x) + 0{x) , a{x) = Aq — -a^x + 0{x^) . (46) 

o 



In the IR, we use the results in equation (43) (up to 0{u^)) as our shooting conditions (we have the freedom 
of choosing the parameter 61, which is always less than zero). We can then extract the values of aq and oi 
at the AdS boundary for our solution. Note that the choice of Bq will determine the choice of ^o- We 
can choose to eliminate the constant A^ by appropriately rescaling our coordinate z. In particular, if we 
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solve our equations with Bq = 0, we can extract the constant £^nd then simply perform the following 
transformation to eliminate it from our metric: 



(47) 



Finally, as indicated in ref. [H], there are constants of the motion regardless of the choice of bi (and 
subsequently, the choice of (ao,ai)), given by: 

-\/6ai 



0.1493 , 



2 In(ao) 



'1.4696 . 



(48) 



These constants, for a given flow, set the size and position of the domain wall separating the two AdSs 
regions of the geometry. 

It is instructive to look at the results to get a sense of what a real domain wall looks like in these 
solutions. We plot (in the original radial coordinate r) the metric function A{r)/r, and also — y4(r)' in 
figure [2] The wall noticeably gets thinner as it moves toward the UV. While it would be instructive to 



A{r) 



d A(r) 





(b) 



Figure 2: Samples of the domain wall behaviour. The lowest (blue solid) curve corresponds to b\ — >■ —3), next lowest 
(green dot-dashed) curve is with bi = —2, and the top (red dashed curve) is with b\ = —1. The UV is to the right, 
and the mass increases with increasing (toward the positive) b\. 



track the wall toward the UV for very large values of the mass, it becomes harder to control the numerical 
accuracy in that regime. 

3.3 The Holographic Entanglement Entropy 
3.3.1 The Box 

The metric we need is given by: 

dsl^ = :^e2(-(2)--(0)) {-de + dxl + dxl + dxl) + ^dz^ . (49) 
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It is important to notice that we are working here in a rescaled z coordinate, z is the coordinate of the 
previous section, with z = z exp(a(0))/i?. This will mean that the mass of the chiral multiplet (controlled 
by the scalar see equation (41 1 ) will now not be set by oq, but by e^'^^ag. 



We consider a strip of infinite extent in the (2:2,2^3) directions and finite extent in the x\ direction. 
We take as coordinates and ansatz for our surface embedding: 

i\=x\ , S,2=X2 , 6 = 2:3 , z = z{xi) . (50) 
The area of the surface is then given by: 

re/2 p3(a(2)-a(0)) , . 1/2 

Area = i?3L2/ dx,- . f 1 + e-2(a(£)-a(o))^,(^^)2y _ (5^) 

J-i/2 Z-i \ ) 

Since there is no explicit dependence on the coordinate xi, a constant of the "motion" is: 

'dz'[x^) z(a;i) Vl + e-2(«(5)-<»(o))z'(a;i)2) ' 

We can evaluate this constant at the turning point x\—Q (with corresponding value z(0) = z*) and rewrite 
our area integral in terms of z: 

rz, 2(a(i!)-a(0)) 

Area = 2R^L^ \ dz , (53) 

J, ^3 /l_e-6(a(i;)-a.)4 



where we have defined our UV cutoff e as z(±£/2) = e. Given the solution for the background, this expression 
for the area can be simply integrated numerically. Furthermore, we can write in terms of t as: 

= / / ■ (54) 

For the case of pure AdS (which simply corresponds to taking a{z) ^ ~ a(0) in our equations), equa- 



tion (53) has a single divergent term near the AdS boundary proportional to e^^. For the fiow geometry, 
the non-trivial behavior of a{z) produces a new divergence in addition to the pure AdS one. Expanding the 
exponential in the numerator, we find that near the AdS boundary, the divergent terms produced by the 



Areauv = 2R'l' I = R'l' ( \ + '^e^'^^'^ In f |, ) ) ■ (55) 



integral are: 



The ln(e) divergence is new. Such new divergence terms associated with the function a{z) are expected in 
rf > 4. 

We will only be interested in the finite contribution (beyond this UV divergent term) to the entan- 
glement entropy, and we will denote this finite contribution by s such that: 
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-0.04- /// 




Figure 3: Blue solid curve is pure AdSs result (corresponds to 61 
and the red dashed curve is with 61 = —1 



—00), green dot-dashed curve is with 61 — —2, 



To perform these integrals numerically, we find that using e = 10"''' and using the coordinate 
y = gives reliable results. We present some results in figure [s] The key result is that wc find that the 
entanglement entropy for the flow geometry asymptotes (for large i/ R) to a constant value not equal to zero, 
as was predicted by the sharp domain wall analysis. The value of the asymptotic constant, which we denote 
by Sooi can be understood as stemming from the parts of the surface in the interpolating region from one 



AdS geometry to the other. With our normalization of s as defined in equation (56), the interpolating part 
of the surface contributes as: 



dz- 



„2(a(£)-ao) 



(57) 



Taking cxd, we find that the finite contribution of this area matches extremely well the finite value 

for the entanglement entropy curves in figure [3] We present the dependence of Soo on the multiplet mass 
in figure |4] We see that for large mass, the value increases as predicted, while for small mass, it actually 
decreases. We can understand this behavior as follows. The sharp domain wall analysis of the previous 
section applies when the mass is large and correspondingly the domain wall's thickness is small. For small 
mass, the wall is fat, and so we can expect deviations (in this case, very large) from our predictions. As the 
mass increases, the sharp domain wall approximation is approached, and our results fit better. Although we 



0.20 
0.15 
0.10 
0.05 
0.00 




„a(0) 

0.8 1.0 1.2^ " 



Figure 4: Value of the asymptotic {i/R — > 00) entanglement entropy for changing multiplet mass. Recall that the 
multiplet mass is proportional to e^'^'ao. 
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have a clear understanding of this asymptotic behavior from our sharp domain wall analysis, where we argued 
it mixes the UV and IR physics around the domain wall, in this example we can completely characterize the 
behaviour in terms of the chiral multiplet mass. We define a function as follows: 

f] = l]o(e"^°'ao)"' , (58) 

where f2o ~ 0.6592i? and is a constant which we found from fitting to the data. Using this function, we 
write the entanglement entropy as: 

AGfs = Area = + ^_^^^Mo) + 2RL\s^,-^ , (59) 

where Sadj now asymptotes to zero. We show this in figure [Sj A nice feature of the function is that it 
captures the IR/UV domain wall physics in terms of a single length scale which naturally normalizes the 
new UV divergence. This length scale emerges naturally from our holographic calculation, and it would 
be difficult (yet worthwhile) to independently calculate it in the dual field theory. Furthermore the non- 
UV divergent parts of the entropy now asymptotes to zero for large £, as appropriate for the entanglement 
entropy for the box region in this limit since there will only be the region A remaining. 




Figure 5: The adjusted entanglement entropy. Blue solid curve is pure AdSs result (corresponds to bi — oo), green 
dot-dashed curve is with &i = —2, and the red dashed curve is with bi = —1. The curves do not perfectly overlap 
when inspected at a higher magnification. 



3.3.2 The Ball 

The metric we need is given by: 

We consider a three dimensional ball with radius i, so take as coordinates and ansatz for our surface 
embedding: 

(i=r, 6 = ^ , 6 = </> , z = z(r) . (61) 
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The area of the surface is then given by: 



ft „2 3(a(i)-a(0)) , . 1/2 

Area = AnR" / dr ^ (l + e'^i-H^-'^W) z' {rf) . (62) 

/g Z[T) \ / 



Since there is an exphcit r dependence, there is no conserved quantity (as was the case for the strip), so the 
extremal area must be found by solving the equations of motion derived from the above action: 

2g-2(a(i)-a(0)) / 3 



z"(r) + z'{rf +{^- 8e2«(")z(r)a'(z) z'(r) 

r \z{r) J 



2 



with the boundary conditions: 



2 3e^(°(^)-°W)(l-2e^"(")^(r)V(g)) _ 

+ z [r) + — J (,Doj 

r z(r) 



z(0) = , z(^) = e . (64) 



In order to study the UV divergences of equation (62), it is convenient to define a new coordinate: 

y{r) = j^/P + e'^-r'^ . (65) 

The reason for this particular choice is that, in the pure AdS geometry, the embedding z{r) has a solution 
given by z(r) — £y{r). In terms of y, the embedding has asymptotic solution given by: 

ziy) ^ey + z^y' - ^^^"^^^^y \n{y) + ■■■ , (66) 

where Z2 is a constant chosen such that the solution satisfies z'{y = 1) = 0. 
Let us define as before a quantity s as follows: 

Our numerical results for s are shown in figure [6] As expected from our sharp domain wall analysis, the 
asymptotic behavior is dominated by an i'^ behavior. Furthermore, upon subtraction of this behavior, we 
find our expected In(^) behavior (see figure [7|. In fact, the large i behavior of the adjusted entanglement 
entropy is very well approximated by the expression: 



hms«a(-j +61n(-)+c. (68) 



Comparing to equation (18), the coefficients a and c should be interpreted as being a complicated function 
of the domain wall data (and hence the mass term), whereas the coefficient b should simply be related to the 
central charge of the IR theory. We show the dependence of a, b, and c on the mass in figure [8j First, we 
see that the coefficient of the area term a decreases with increasing mass, which deviates from our domain 
wall analysis. However, our mass range displayed here is actually in the same region where the analogous 
results for the box were decreasing. This is all traceable (as it was then) to being far from the thin domain 
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Figure 6: Results for the ball with radius I. Blue solid curve is for —bo = 3, green dot-dashed curve is for —bo = 2, 
and red dashed curve is for —bo = 1. 

P 

s - a — 




-0.10^ 



Figure 7: Entanglement entropy without the leading contribution. Blue solid curve is for —bo = 3, green dot- 
dashed curve is for — 6o = 2, and red dashed curve is for —bo = 1. 



wall regime. In fact it is apparent that our results are just about to turn around, as happened for the box. 
We may ask whether we can capture the behavior of this coefficient with a function Q, defined in a similar 



way to that which we did for the box (see equation (58) and discussion below it). This turns out to be the 
case, and we compare the function il for the strip and for the ball in figure [9] ( where we again let the data 
determine the value flo) and we find that the results are remarkably close, suggesting that indeed this is a 
physical length scale in the theory. Since the only other scale in the theory is the multiplct mass, we expect 
that il is set by the inverse of the value of this mass after it has renormalized under RG flow to its value set 
by the scale at the domain wall. 

Finally, it is interesting to check whether the coefficient of the \n{£) term of the entanglement 
entropy indeed matches our expectation that for large enough £ it should be equal to —R'l^/2R^ (compare 



the definition of s in equation (67), and b in equation (68) and (18)). We present this check in figure 10 
Indeed, we find that our results are very close to the expected value, and in fact the results get better with 
larger mass. It is not entirely unexpected that the results deviate somewhat for smaller mass. With smaller 
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Figure 8: Values of the fit coefficients of equation (68 1 



masses, one needs to go to larger £ (and larger AdS radial position) to enter the IR AdS, and larger £ is 
increasingly hard to explore numerically. 

4 The Three Dimensional Holographic RG Flow 

4.1 The Holographic Dual Gravity Background 

We consider the 2-fl flow of ref. 525|. It is the analogue of the Leigh-Strassler flow we studied earlier. The 
IR fixed point is the Af — 2 fixed point studied in [321 ED] ■ Starting with an ansatz for the metric of the 
form: 

dsls = e^^^'') {-dt^ + dxl -f dxl) + dr^ , (69) 
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Figure 9: Comparing the prediction for the function f2 for the strip and the ball. Blue crosses are the ball results 
and the red squares are the strip results. 
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Figure 10: The coefficient of ln(£) for large £. The solid line is the value of —Rfj^/R^. 



the flow equations are given bj{^ 

A'ir) 
p'{r) 

X'ir) 



1 (cosh(2x) + 1) + (cosh(2x) - 3) 
16^ p 
1 (p« - 3) sinh(2x) 



4i? p2 

To study the UV asymptotics, we define a coordinate z via: 



(70) 



(71) 



^We use a different definition of the AdS radius R than the authors of ref. |2S] . The reason for this is such that the asymptotic 
behavior in terms of the coordinate r are the same for our previous analysis. 
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and find that the fields have the following asymptotics: 

A{i) - - ln(z)-^ (12a? + a2)+ 0(^3) ^ 
x{z) — za^ — 6aiao2^ + 0(z"^) , 

a{z) = zai + (^al - z"^ + 0{z^) . (72) 
In the IR, we define a coordinate u via: 

u = e^'-/^ , (73) 
where A = 3'^/'' (vTY— l) /4. The asymptotics of the field is given by: 

A{u) = ^ ^ In(fi) + So + 0{u^) , 



Xiu) = ^cosh \2) + bou + 0{u^) , 



a{u) = Un{i)+'^^^^b^u + 0{u^) . (74) 
Note that this gives that i?iR = 2^/3^/". 
4.2 The Entanglement Entropy 

We follow a similar analysis as before for calculating the entanglement entropy and in light of our detailed 
calculations presented for AdSs, we will be brief here. For the strip, we define a quantity s via: 

/ f?2 r \ 

4GnS ^ Area = + 2RLs , (75) 



and we plot the result in figure 11(a) Similarly for the disc (of radius i), we define s via: 

4:GnS = Area = 2nR'^ + s] , (76) 



and we plot the result in figure [lT(b) As expected, there is no new divergent term, and the entanglement 
entropy has a non-zero coefficient associated with the area of the boundary for large £. Furthermore, fitting 
the large i behavior for the disc result to the functional form a£ + b gives a result for b of approximately —0.79 
which is close to the value of — -Rjr, the value we would have predicted in the IR AdS from the thin domain 



wall analysis (see equation (21 1 and discussion below it). Using the strip data, we calculate the dependence 
of the coefficient on the mass and present it in figure [T2j We find that the dependence is linear in the mass 
range explored. The linear behavior again suggests that we can interpret the coefficient in terms of the mass 
of the multiplet. The slope is negative. We expect that we are somewhat away from the mass regime where 
we connect with the thin domain wall analysis for physics controlled by the domain wall region. The regime 



of larger mass ought to show that the coefficient increases towards positivity (see equations (22) and (30)). 
It is currently too difficult to extract reliable numerical results in that regime however, and so we cannot 
test this as we did for the case of one dimension higher. 
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Figure 11: Results for the strip of width £ and the disc of radius I. Solid blue is the pure AdS4 result, green dot 
dashed is in a background with — 6o = 0.01, and red dashed is in a background with —bo = 0.001. 
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Figure 12: The coefficient of i for large £. The blue circles are the numerical results, while the red dashed curve is 
the linear best fit. 

5 Conclusions 



We were able to uncover how the holographic entanglement entropy formula encodes a number of field theory 
features along an RG flow using the key feature of an holographic RG flow, its domain wall. In the thin 
wall limit where we study an idealized flow between two AdS regions (and hence a flow between fixed point 
theories) we were able to extract, in various dimensions, rather pleasing formulae showing how the field 
theory data appear separated out according to the various natural scales in the problem. Using two known 
RG flows between fixed points, we were able to test our formulae in real examples and found that where 
we could, several of the features we expected to be robust were confirmed. This included the form of the 
growth of the entropy upon approach to the IR, which we found defined a new area law supplementing the 
ones known for the UV behaviour. We were also able to see (in the thin domain wall limit) how the entropy 
naturally encodes (in its asymptotics) the (renormalized) length scale corresponding to the presence of the 
relevant operator that was switched on. 
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Our overall goal is to understand the entanglement entropy in this holographic setting well enough 
to use it as another diagnostic tool for studying the properties of strongly coupled field theories, and so 
we consider this study of RG flow to be a useful step along the way. Having analytically characterized, as 
we have done, the kinds of behaviour that can appear (in various dimensions), we expect that holographic 
studies of the entanglement entropy in more complicated theories will be aided by our results. Many such 
examples will, by their very nature, be only accessible with numerical approaches, and so analyzing of the 
results for the entropy may be subtle. This is where we expect a lot of the intuition here to help bring 
things sharply into focus, since such features as the entropy's approach to the end of the flow have now been 
unpacked in terms of the general features of holographic flow geometries. 
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